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ABSTRACT 

Aims. We investigate the effect of the physical environment on water and ammonia abundances across the S140 photodissociation 
region (PDR) with an embedded outflow. 

Methods. We have used the Odin satellite to obtain strip maps of the ground-state rotational transitions of ortho-water and ortho- 
ammonia, as well as CO(5-4) and 13 CO(5— 4) across the PDR, and H^ 8 in the central position. A physi-chemical inhomogeneous 
PDR model was used to compute the temperature and abundance distributions for water, ammonia and CO. A multi-zone escape 
probability method then calculated the level populations and intensity distributions. These results are compared to a homogeneous 
model computed with an enhanced version of the RADEX code. 

Results. H 2 0, NH 3 and 13 CO show emission from an extended PDR with a narrow line width of ~3 km s _1 . Like CO, the water line 
profile is dominated by outflow emission, however, mainly in the red wing. Even though CO shows strong self-absorption, no signs 
of self-absorption are seen in the water line. H' s O is not detected. The PDR model suggests that the water emission mainly arises 
from the surfaces of optically thick, high density clumps with HflD^lO 6 cirT 3 and a clump water abundance, with respect to H2, of 
5x 10~ 8 . The mean water abundance in the PDR is 5x 10~ 5 , and between ~2x 10~ 8 - 2x 10~ 7 in the outflow derived from a simple two- 
level approximation. The RADEX model points to a somewhat higher average PDR water abundance of lxl0~ 8 . At low temperatures 
deep in the cloud the water emission is weaker, likely due to adsorption onto dust grains, while ammonia is still abundant. Ammonia 
is also observed in the extended clumpy PDR, likely from the same high density and warm clumps as water. The average ammonia 
abundance is about the same as for water: 4xl0~ 9 and 8xl0~ 9 given by the PDR model and RADEX, respectively. The differences 
between the models most likely arise due to uncertainties in density, beam-filling and volume filling of clumps. The similarity of 
water and ammonia PDR emission is also seen in the almost identical line profiles observed close to the bright rim. Around the 
central position, ammonia also shows some outflow emission although weaker than water in the red wing. Predictions of the H2O 
1 1 ~ lo.i and I11 -Oo.o antenna temperatures across the PDR are estimated with our PDR model for the forthcoming observations 
with the Herschel Space Observatory. 

Key words. ISM: abundances - ISM: individual (Sh 2-140) - ISM: molecules - Submillimeter - Line: profiles - Line: formation 

1. Introduction I2003I) have both shown very low a bundances of molecular 

^ , , , . . , . „ ,. oxygen, 5xl0~ 8 dLarsson et a l„ 2007), and water, 10~ 8 -10~ 9 

The chemistry and the evolution of the interstellar ^medium de- (Qlo fsson et all \m0$,]n the extended cold clouds in contrast 

pend to a high degree on the available amount of carbon and tQ ^ dictions from chemical models However, at smaller 

oxygen. The carbon-based chemistry is fairly well understood, scales ^ wam Qr shocked ionfj md outnows ^ where the den- 

but the oxygen chemistry is not. Two molecules vital to the tin- ^. M hej . than 1q5 _ 1Q 6 cm -3 SWAS ( MelnickEH 

derstanding of the oxygen chemistry are water and molecular 2000a) and Qdm ( Ipersson et ail 50O7h have also shown lo- 

oxygen, since some chemical models predict that a large frac- cM increased water abundance by severa i orders of magnitude, 

tion of oxygen could be in the form of these molecules. 1Q - 6 _ 10 _ 4 with ct tQ molecular hydrogeri . The analysis of 

• ? f^Tl § rou "?- stat A e rotational water transition ^ ^ ^ somewhat difficuh and uncerta i n due to 
with NASA s Submillimeter Wave Astronomy Sate llite (SWAS ; observations of a si le line and k beam-widths, but is in 
i Melnicketalj |2000bJ) a nd the Odin satellite (|Nordhetal| agree ment with observations of a wealth of highly excited wa- 

Send offprint requests to: carina.persson@chalmers.se ter vapor transitions with ISO (the Infrared Space Observatory, 

* Odin is a Swedish-led satellite project funded jointly by the |Cernicharo & Crovisieij |2005|). These locally high abundances 
Swedish National Space Board (SNSB), the Canadian Space Agency support the predictions that water should be an important cool- 
(CS A), the National Technology Agency of Finland (Tekes) and Centre ing agent in the star-formin g process, especially at high tem- 
National d'Etudes Spatiales (CNES). The Swedish Space Corporation peratu res and densities (e.g. iNeufeld et all 1 1 995t iBergin et all 
was the prime contractor and also is responsible for the satellite opera- |2003l) . The strong dependence of the water chemistry on the tem- 
tion. 
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perature changes seen in star-forming regions also make water 
an excellent tracer of shocks and outflo ws from protostar s still 
hidden within their dusty envelopes (e.g. lBergin et al.Lll998l) . 

The intention of this paper is to continue the previous work 
with ISO, SWAS, and Odin, and also to be a valuable input to 
future observations with the Heterodyne Instrument for the Far 
Infrared (HIFI) on board the HerscheQ satellite to be launched 
in 2009. Our aim is to investigate how the water abundance de- 
pends on the environment in a star-forming region. The high 
spectral resolution (0.5kms~') and the smaller beam-size of 
Odin compared to SWAS (2.'1 vs. 3 .'3x4^5), which enable higher 
spatial resolution maps, allows a more detailed analysis of the 
line profile observed at each map position. We can also obtain 
an estimate of the beam-filling by comparing SWAS and Odin 
line profiles if the same object has been observed. If not prop- 
erly taken into account the beam-filling may otherwise influence 
the results by orders of magnitude. 

In addition to the water observations a second goal of our 
work is to perform and analyse observations of the important 
ground-state rotational transition of ammonia. While ammonia 
has been extensivel y observed in the inversion transitions a t cm- 
wavelengths (e.g. lUngerechts et al.iri986tlZhou et al.lll993l) . the 
lo-Oo transition, however, cannot be observed from ground 
due to the opaque terrestrial atmosphere. The first observa- 
tion of this transition was perf ormed by the Kuiper Airborne 
Observatory ( Keene et aTlll983l), and was later followed up with 



20031: iLiseau et all 120031; blofsson et al.1 



Odin dLarsson et al 
l2007h . 

Our chosen target is the well known molecular cloud 
and bright nebula Sh 2-140, more commonly known as S140 



dSharplessl 19591 iBlair et aU 19781). It is located at a distance 
of 910 pc dCrampton & Fished 1 1974 in the large L1202/L1204 
dark cloud. This region is illuminated from the south-western 
side by the nearby B0.5V star HD 21 1880, creating a bright 
HII region and an almost edge-on PDR adjacent to the molec- 
ular cloud and extending deeply into the cloud. The core in 
S140 also contains a small cluster of deeply embedded young 
stellar objects which illuminate the cloud from within. One 
or more of these early B-stars are undergoing some mass-loss 
process and molecula r outflows are observed from this dense 
and warm reg i on (e .g. Rouan et al. 1 1 1 9771 lHavashi et al J. fl987t 
iMinchin et aU 11993 ;lPreibisch & Smithl |2002|) . The brightest 
infrared source in S140, called IRS1, is on small scales com- 
posed of several sources, two of which might be powering 
bipolar molecular ou tflows in almost perpendicular directions 
dTrinidad et al.L 120071) . The increased temperature and density 
in both the PDR and the outflows are likely expected to produce 
a higher than average water abundance. For these reasons, SI 40 

is a promising target for our purposes. 

Both ISO dBoonman & van Dishoecld l2003h and SWAS 
dSnell et all lToOOa; Ashb v et al.1 , 2000a b; Franklin et all 120081) 



have observed water in S140 and several at tempts have 
been ma de to model th e se obs e rvations (e.g. lAshbv et all 
2000ab; Boonman et a 1 120031; iPoelman & Spaansl 120051 
20061) . However, they all have had difficulties to explain the un- 
usual, narrow li,o-lo,i emission line with no or little evidence 
of self-absorption. In order to obtain more information about 
the water emission we have conducted strip map observations 
of S140, starting deep within the molecular cloud and reaching 
outside the bright rim (see Fig. [T]). Strip maps are also simul- 
taneously performed of the previously unobserved NH3 (lo-Oo) 
transition. 




http://herschel.esac.esa.int/ 



Fig.l. A CI map (4 92 GHz) of S140 observed with JCMT 
( Min chin etail Il993b . The Odin 2'A beam at the five strip po- 
sitions is also shown. Position 1 in NE to position 5 in SW. 



To aid the analysis of water and ammonia, we have in addi- 
tion performed Hj 8 observations in the central position, strip 
maps of CO (5-4) and 13 CO (5-4) with the Odin satellite, and 
a 13 CO (1 -0) map with the Onsala 20-m telescope. 



2. Observations 

All transitions, except 13 CO(l -0), were acquired with the Odin 
satellite from 2003 to 2006 with a total of 561 orbits (see 
Table [T). The calibr ation procedure (t he chopper wheel method) 
is described in Olberg et al.J d2003l) . The Odin 1.1m offset 
Gregorian telescope has a circular beam at 5 57 GHz with a Ful l 
Width of Half Maximum (FWHM) of 2.'1 dFrisk et all |2003). 
Being outside the atmosphere, and with an exceptionally high 
main beam efficiency, rj m \, = 0.9, our intensity calibration is very 
accurate. The intensity scale in the figures is expressed in terms 
of antenna temperature Tl. In all calculations, however, the 
main beam efficiency is properly taken into account. The recon- 
structed pointing offset was <15" during most of the time. The 
simultaneous observations of H2O and NH3, and of 13 CO, CO 
and H2O guarantee their pointing to be identical. 

Three different tunable submm receivers were used, hav- 
ing average single-sideband (SSB) system temperatures of about 
3 300 K. They were used in combination with a hybrid autocorre- 
lator spectrometer (300 -400 MHz working bandwidth) and an 
acousto-optical spectrometer (1 GHz working bandwidth), both 
with a channel spacing of ~0.3kms _1 . All observations were 
performed in the sky-switching mode with reference sky beams 
at 42° distance and 4.4° FWHM. To correct for ripple, observa- 
tions of a reference off position 30' east of S 140 were conducted 
as well. 

The five point strip maps have 60" step size with strip centre 
coordinates R.A. 22 1 M7M2 s .O, Dec. +63°03:45'/0 (B1950) at the 
cloud core (our position 3, Fig.Q]). Position 1 is in NE and posi- 
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Table 1. Observed molecular transitions. All lines are observed 
with the Odin satellite except the 13 CO(l -0) map which is ob- 
served with the Onsala 20-m telescope. 



Species 


Transition 


Freq. 


K 


n", 

cnt 


Int. time 






[MHz] 


[K] 


[cm ] 


M 


H 2 


1 1,0 _ lo, 1 


556 936.0 


61 


3xl0 8c 


41 


H> 8 


1 i,o — lo, 1 


547 676.4 


60 


3xl0 8c 


91 


CO 


5-4 


576267.9 


83 


2xl0 5c 


9 


13 co 


5-4 


550 926.3 


79 


2xl0 5 


9 


13 co 


1-0 


110 210.4 


5 


7xl0 2 


2 min/pos 


NH 3 


lo-Oo 


572498.1 


29 


5xl0 7c 


16 



"Upper state energy. 

''Critical density calculated for a temperature of 40 K. 

' Can be considerably reduced because of radiative trapping. 



tion 5 in SW. The frequency scale is set relative to a source LSR 
velocity of -7.5 km s . 

In addition, the Onsala 20-m telescope was used in 2006 
for a 13 CO(1-0) 10x11 point map of the S140 molecular 
cloud with a spacing of 30". The SSB system temperature was 
~500K, the main-beam efficiency = 0.5, and the FWHM beam 
at 110.201 GHz is 34". 



3. Results 

The observed Odin spectra of CO, H2O, and NH3 are shown 
together with 13 CO(5-4) in Figs. [2]— 14] to allow a compari- 
son of the line profiles of the different molecules at each posi- 
tion. The water and ammonia spectra are also shown together in 
Fig. [5] No H^O emission was detected in the central position 
(Fig. IB.ll on-line material) at an rms noise level of 8 mK. Tables 
with lcr rms, integrated intensities, peak antenna temperatures, 
line widths and amplitudes of Gaussian fits, and centre veloc- 
ities for species observed with Odin are found in the on-line 
Tables IB. 1 1 — IB.41 The on-line material also includes figures of 
Gaussian fits to 13 CO(l -0) and 7 = 5 -4, H 2 and NH 3 (li -0 ) 
(Figs.lB~4l-lBJl), 

Below we will discuss the observations and some straight 
forward results, while the model results for the narrow PDR 
component are presented in Sect.|4]and|5] 

3.1. Carbon monoxide 

The CO(5-4) and 13 CO(5-4) observations show a narrow line 
from the PDR in all positions superimposed on a broader outflow 
feature seen in CO at positions 1-4 and in 13 CO at the centre 
position (Fig. [2]). The emission peaks about 20" from the central 
position towards position 2 deeper in the cloud. A comparison of 
the 12 CO and 13 CO line shapes suggests strong CO(5-4) self- 
absorption by lower excitation foreg rou nd gas. This is also dis- 
cussed by e.g. M inchin etail d 19931) and iBoonman et alj d2003l) 
where the higher /-level CO lines are strongly self-absorbed. No 
strong self-absorption is visible in 13 CO. 

The mean opacity of 12 CO and 13 CO can be determined by 
comparison of their antenna peak temperatures if the isotopo- 
logue abundance ratio R is known. Assuming that the opacity of 



12 CO is larger than the opacity of 13 CO by a factor R, we will 
have 

2Xi3 = 7(0(1 - e - T ")77bf,i3 = 7(0(1 -e- T ")rj bm 
T ln J(T£) (1 - e-n) 77bf,i2 J{Tlb (1 - 77bf,i2 ' 

where the radiation temperature, J(T ex ), is a function of the ex- 
citation temperature, T ex , 

The beam-filling factor rj^f = 2 /(9 2 + 9\), assuming that 
both the source brightness distribution and the antenna response 
are circularly symmetric and Gaussian. 

Assuming that the excitation temperatures are the same for 
both species, co-spatial emissions and using R = 70, the opaci- 
ties of the narrow component of 12 CO and 13 CO are estimated to 
be < 70 and <1, while the opacities in the broad outflow compo- 
nent are <7 and <0.1, respectively. Due to the self-absorption 
of 12 CO the opacities are upper limits. 

The size of the emission region can be estimated by compar- 
ing Odin and S WAS antenna temperatures 

r l,od,n = (33 x 60 x 4:5 x 60) + 6 2 

T hwAS (2.'lx6O) 2 + 2 ' {y> 

where 8 is the effective circular Gaussian source size in arcsec- 
onds. For the central position the ratio of the Odin and SWAS 
13 CO peak antenna temperatures is ~2.5 which is equivalent to 
a mean source size of ~100". 

A second approach to estimate the source size is performed 
by plotting the antenna temperatures in our Odin strip map vs. 
position, and fitting a Gaussian. After deconvolution with the 
126" Odin beam, the resulting one dimensional source size is 
about 106". Together with the mean source size of ~100", this 
implies a relatively circular source. 

The excitation temperature can be estimated by plotting a 
rotation diagram of several transitions. The level populations are 
given by the Boltzmann equation and the excitation temperatures 
for all the energy levels are assumed to be the same. The beam- 
filling corrected integrated intensities of all lines are plotted as a 
function of t he upper state energ y in a semi-log plot based upon 
the relation (Perss on et al. , 2007) 

in -i^L fr b d, = ln^I-^. (4) 
hc^guAuJ Q(T) kT ex 

where k is the Boltzmann constant, v„/ is the frequency of the 
transition, h is the Planck constant, c is the speed of light, A„/ is 
the Einstein A-coefficient for the transition, 7\, is the brightness 
temperature (observed antenna temperature corrected for beam- 
efficiency and beam-filling), g u and E u are the statistical weight 
and energy of the upper state, respectively, Q(T) is the partition 
function, T ex is the excitation temperature for the transition, and 
A^rot is the total column density of the species obtained from 
the rotational diagram. As customary the frequency axis v has 
been converted to a velocity axis v using the speed of light. A 
least-squares fit to the data will produce a straight line with slope 
-I/Trot- If we extrapolate the line to E„ = K, the total column 
density is found from the intersection of the y-axis. 

Our rotation diagram of the broad outflow 13 CO compo- 
nent in the central position (on-line Fig. |B.2t , including the 
Onsala 20 -m 7=1-0 transition , the 7 = 2-1 and 3-2 transi- 
tions from Minchi rTet al.l (Il993l) . and the Odin 7 = 5-4 transi- 
tion, produces r R0T = 24+2 K and N RO t = (2.5±0.4)xl0 16 cirT 2 . 
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Fig. 2. Comparison of the CO(5-4) (thick magenta line) and 
13 CO(5-4) (thin black line) line profiles in a strip map from 
NE to SW across the PDR. 



Velocity [km s ] 

Fig. 3. Comparison of the H2O li,o-lo,i (thick magenta line) 
and 13 CO(5 - 4) (thin black line) line profiles in a strip map from 
NE to SW across the PDR. 



All antenna temperatures are corrected with respective beam- 
filling and beam-efficiencies. This temperature is in agreement 
with the previously estimated lower limit of the outflow e xcita- 
tion temperature of about 15-20K dMinchin et all Il993b . The 
low temperature of the outflow indicates that the emission arises 
from g as far behind the s hock front where the gas has already 
cooled dSnell et al.Lll984bl) . 

The excitation temperature of the narrow I3 CO PDR emis- 
sion in the central position is also estimated by means of a rota- 
tion diagram (on-line Fig. IB. 31 ). including the above mentioned 
transit ions with addition of the 7 = 6-5 line from iGraf et al] 
(1993), but excluding 7=1-0 which seems to originate only 
in the colder ga s component. T his produces Trot = 69±27 K, in 
agreement with Mi nchin et al. I (fl993l) who find a PDR excitation 
temperature of about 70 K at IRS and a steep rise in temperature 
to about 250 K at the dissociation front 60" -80" SW of IRS1. 
Our PDR column density, V ROT ( 13 CO) = (3.2+1.8)xl0 16 cirT 2 , 
is relatively uncertain but in agreement with the RADEX r esult i n 
Sect.H V( 13 CO) = 2.8xl0 16 cm 2 , and with lPlume etai](ll994l) . 
3.4xl0 16 ctrT 2 . 



3.2. Water 

The ortho-water 1 1,0 — lo.i ground-state rotational transition is 
detected in all positions except outside the bright rim (po- 



sition 5). Similar to CO and 13 CO, a narrow component 
(Av ~3kms~ I ) superimposed on a broader outflow component 
(Av ~8.5kms ) is observed with emission originating around 
the central position. Like CO, the water line profile is dominated 
by outflow emission, clearly seen in a fit of two Gaussian compo- 
nents in the on-line Fig. IB. 61 The integrated intensity of outflow 
emission is only about 20% for 13 CO compared to 60% for water 
which in addition mainly shows a red wing (Fig.[3j. This appar- 
ent strong outflow emission may not, however, reflect the true 
relative abundances of the PDR and outflow components if the 
line is optically thick. The opacity is several times higher in the 
PDR than in the outflow and, depending on the opacity, this may 
lower the fraction of outflow emission. The I3 CO emission is on 
the other hand almost optically thin which gives more accurate 
relative amounts of the respective emissions. 

The non-detection of Hi 8 O can be used to set upper limits on 
the opacity of the main water line using Eq.Q] However, the cal- 
culation is not as straightforward as for CO. Because of the high 
critical density (the ratio of the spontaneous de-excitation and 
collision coefficients) of the H2O 557 GHz line, the emission of 
this transition is highly sub-thermal and has a high opacity even 
for rather low abundances. This high optical depth allows the 
excitation temperature to be enhanced by radiative trapping. As 
a result, the corresponding lines of H2O and H^O have quite 
different values of excitation temperature. Their excitation tem- 
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Fig. 4. Comparison of the NH3 lo-Oo (thick magenta line) and Fig. 5. Comparison of the H2O (thick magenta line) li,o~ 
13 CO(5 -4) (thin black line) line profiles in a strip map from NE 1q,i and the NH3 (thin black line) Iq-Oq line profiles in a strip 



to SW across the PDR. 



map from NE to SW across the PDR. 



peratures will also be much lower than the kinetic temperature. 
Using the RADEX code in Sect. [4] for a kinetic temperature of 
55 K and a molecular hydrogen density n(H2) = 4x 10 5 crrT 3 , we 
find excitation temperatures of 8.2 and 6.3 K for the 1 1,0 - ln,i 
transitions of H2O and H^O, respectively. If the abundance ra- 
tio [ 16 0/ 18 0] = 330 in water dPersson et aUl2007l) and if the two 
isotopologues are distributed similarly, then the 5cr upper limit 
on the intensity of the Hi 8 O line implies that the optical depth of 
the H2O line is < 60. If the assumed density is lower by a factor 
of 10, the difference in excitation temperatures will increase and 
change the upper limit to <120. 

In addition, the high H2O opacity often produces self- 
absorbed spectra. Our observations, however, do not show any 
obvious signs of self-absorption, in contrast to our CO lines. 
When comparing the centre velocities of the narrow components 
of water and CO to 13 CO, we find that the velocities of water 
and 13 CO are well matched (Fig. [3] and Table IB~2l-lB~3l, while 
CO and 13 CO in positions 1-4 are not (Fig.|2]i. This supports 
our conclusion that there is no severe self-absorption in H2O in 
contrast to observations of other sources in which corresponding 
self-abso rptions were seen in C O(7-6) and in ISO measurements 
of H 2 dBoonman et aUl2003l) . 

The lack of self-absorpti on points to a rathe r constant exci- 
tation across the region (cf. ISnell et all Il984ah . It also implies 
a lower water abundance in the blue than in the red outflow. A 



higher water abundance in the red wing may be caused by out- 
flowing gas pushing into the ambient high density cloud thereby 
producing water. The blue outflow, on the other hand, is leaving 
the molecular cloud, and expands into an ionized region with low 
density gas with no or very low w ater production as a r esult. This 
is in agreement with the results in Fra nklin et al.l d2008l) who in- 
vestigated water abundance in molecular outflows with SWAS 
and find about five times higher abundance in the red outflow 
compared to the blue. 

The sizes of the emitting regions are estimated in the same 
way as for I3 CO (Eq.0. For the central position the ratio of the 
Odin and SWAS peak antenna temperatures is ~2.4, equivalent 
to a mean source size of ~100". This confirms previous sug- 
gestions that the water emission does not fill the SWAS beam 
dBoonman et all [2.0031) . To obtain estimates of the size of each 
component, we use the amplitudes of the fitted Gaussians from 
the Odin and SWAS spectra. The ratios of the PDR and out- 
flow components are ~2.2 and ~2.6, which corresponds to mean 
source sizes of 120" and 85", respectively. 

Gaussian fits to the amplitudes in each position in the Odin 
strip map results in one-dimensional source sizes of >150" and 
<50" for the PDR and outflow, respectively. Taken together, 
these two calculations of source sizes suggest that the water 
emission from the PDR is more extended in the NE - S W direc- 
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tion (same direction as our strip map), while the contrary applies 
to the outflow, indicating a NW - SE elongation. 

A third approach to estimate the beam-filling can be achieved 
by using the radiative transfer equation. With a constant source 
funct ion, the solution of the radiative transport equation is 
(cf. lPersson etaTl 120071) 



T* A = T b rj mh T] bf = (J(T ex ) - J(T conl )) (l-e T ) rj mh T] bf , 



(5) 



where J(T cont ) is the background continuum radiation. 

For an optically thick line and neglecting the background ra- 
diation 



J(T ex )rimb Tbf. 



(6) 



If r ex is known, then the beam-filling can be estimated. Using 
r ex = 8.2 K obtained from RADEX and Equation [6] give a beam- 
filling factor that corresponds to a PDR source size of about 
110", which is close to the source size of 120" obtained from 
the relative temperatures from Odin and SWAS. This implies 
that the gas is somewhat clumped, and that the surface filling 
fraction of the high density PDR gas is ~84%. The clumpy 
structure of the clo u d is confirmed by previous observations 



Snell et all Il984albt IZhou et all Il993t iMinchin et all Il993t 
Havashi & Murata, 1992; Ash bv et aUl2000al) and our model re- 
sults in Sect. [5] The amount of surface clumping is, however, not 
possible to determine very accurately due to the uncertainty of 
the excitation temperature which is affected by the assumed den- 
sity. With T ex = 9K and 10 K the same calculations give source 
sizes of 86" and 69", with corresponding surface filling frac- 
tions of 50% and 33%, respectively. Including the background 
radiation in Eq.|5]also introduces additional uncertainties. 



3.3. Ammonia 

A complex energy level diagram with a wealth of transitions 
covering different temperatures and densities has made ammonia 
a very valuable diagnostic of physical conditions in the interstel- 
lar medium. The lo-Oo ground-state rotational transition, ob- 
served by Odin, has an upper state energy comp arable to the ex- 
tensively observed inver sion transitions (e.g. lUngerechts et all 
1986;|Zhou et al., 1993), but the critical density is about four 
orders of magnitude larger. It is therefore likely that these tran- 
sitions partly probe different gas components. As for water, the 
lo-Oo transition has a high opacity and is highly sub-thermally 
excited in most cases. 

The lo-Oo transition of ortho-NHi is detected in all posi- 
tions except outside the bright rim (position 5). The line profiles 
and centre velocities are very similar to those of 13 CO(5 -4) with 
an exception in position 1 (Fig. |4j». Our comparison of ammonia 
and water in Fig. [5] shows an almost identical line profile at the 
bright rim (position 4), which suggests emission from the same 
gas and velocity fields at this position originating in the PDR. In 
the central position the water emission shows a more pronounced 
outflow in the red wing than ammonia, although very similar in 
the blue outflow. The emission fraction of the outflow is 35% for 
ammonia, while it is almost twice as high for water. This sug- 
gests that the outflow production mechanism is more efficient 
for water than ammonia. The water emission from the outflow 
seems to be almost the same in position 2 and 3, while the ammo- 
nia emission from the outflow decreases in position 2. Deepest in 
the molecular cloud a very weak water line is seen. The ammonia 
emission is, however, considerably stronger than the water emis- 
sion at this position. This confirms prev i ous observations of NH3 
inversion transitions (Zh ou et al l ll993HUngerechts et allfl986h 



who find a NH3 peak further into the cloud at about 1' northeast 
of IRS1. When we plot the peak antenna temperatures vs. posi- 
tion it is obvious that we see additional emission in position 1 . A 
Gaussian fit to positions 2-4 gives a source size of ~130", with 
a peak about 10" into the cloud from the centre. This size is rel- 
atively uncertain due to a fit to only three positions, but agrees 
well with the source sizes estimated from 13 CO and H2O. 

The Odin observations indicate that in addition to the warm, 
dense PDR and a weak outflow around IRS, we also observe 
NH3 from the same cold and extended gas deep in the cloud (po- 
sition 1) as 13 CO(l - 0) at a temperature about 2 - 3 K. This is 
also in agreement with Ung erechts et alj d 1 986b and IZhou et all 
( 1993) who find a NH3 rotation temperature of 40 K in between 
our first and second position. Chemical models predict that NH3 
is not depleted in high density regions (<lxl0 6 cm" 3 ) and is 
therefore relati vely more abundant than other molecu les such as 
water and CO dBergin & Langerl 119971; iBerginl 120031) . 

The surface filling fraction of the narrow PDR component of 
ammonia is also estimated using the radiative transfer equation 
(Eq. |6j. The excitation temperature is approximately the same 
as for water thereby producing similar surface filling factors. 
Together with the similarity of the line profiles, source sizes and 
centre velocities of the H2O, NH3 and 13 CO(5 -4) narrow com- 
ponents, this points to emission from the same high density PDR 
clum ps with the same tem perature of about 70 K for all species 
(cf . iBoonman et all 120031) . 

3.4. Outflow abundances 

The collisional de-excitation of the ground-state rotational tran- 
sition of ortho-water is very low in the density ranges considered 
here due to the high critical density of ~3xl0 8 cm" 3 at a tem- 
perature of 40 K. A simple analytic expression for the antenna 
temperature of a two-level system in the low collision rate limit 
can therefore be used since collisional excitation always results 
in a photon that escapes the cloud. Thus, even thoug h our tran- 
sition may be op t ically t hick, it is effectively thin dLinke et all 
\19T% ISnell etafl. l2000bl) . We u se this approach to estimate the 
water abundance in the outflow (ISnell et all 12000b) 



T b dv = C n Hl ^— N(o - H 2 0)^ exp(-hv/kT K ), (7) 



where C is the collision rate which equals the product of the 
cross-section and velocity. Di viding Eq. j T]) with jV(H2) and scal- 
ing the calculated factor in ISnell et ahl d2000bl) to our derived 
rotational outflow temperature of 24 K (with yT) produces the 
fractional abundance of ortho-water with respect to H2 in the 
outflow 



X(o - H 2 0) = 3.2 x 10 



19 



J 7b dv 
N(H 2 )nn 2 



(8) 



The H2 column density in the outflow is calculated 
from 7Yrot( 13 CO) of the broad component. Using typical 
values of [CO/H 2 ] = lxlO" 4 and [ 12 C/ 13 C]=70 we obtain 
7V(H2) = 1.8 xlO 22 cm" 2 . Together with a range of mean molecu- 
lar hydrogen densities of 5xl0 4 -5xl0 5 cm" 3 we obtain a total 
water outfl ow abundance of ~ 2x 10" 8 - 2x 10" 7 . This is in agree- 
ment with Franklin et al. (2008) who found an outflow abun- 
dance in the blue wing of 2.4xl0" 8 and l.OxlO" 7 in the red wing 
using a kinetic temperature of 30 K and a molecular density of 
10 5 cm" 3 . 

The ammonia abundance in the outflow is estimated using 
the Local Thermal Equilibrium (LTE) approximation. Assuming 
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Table 2. Observations and model calculations for the narrow PDR component in the RADEX homogeneous model." 



Species 


Transition 


Frequency 

V 

[GHz] 


T b 
[K] 


Observations 
Anl 

[ergs s crrT 2 ] 


Av 
[kms" 1 ] 


Ref 


T 

* ex 

[K] 


Model calculations 
t T b Anl 

[K] [ergs s _1 cm" 2 ] 


H 2 


1 1,0 _ lo.i 


556.936 


0.92 




3.1 


c 


8.2 


6.6 


0.92 


6.8 x 10- b 


H' s O 


1 i,o — lo.i 


547.676 


< 0.051 




(3.1)* 


c 


6.3 


0.2 


< 0.05 


< 3.5 x 10" 7 


NH 3 


lo-Oo 


572.498 


1.1 




3.3 


c 


9.0 


U. 1 


[ [ 


9.2 x 10~ 6 




(1,1) 


23.695 


2.9 






d 


54 


0.02 


0.87 


2.5 x 10- 10 




(2,2) 


23.723 


2.0 






d 


47 


0.01 


0.50 


1.7 x 10" 10 




(3,3) 


23.870 


1.4 






d 


-5.1 


-0.25 


2.2 


9.0 x 10" 10 




(4,4) 


24.139 


0.15 






d 


46 


0.001 


0.06 


3.8 x 10"" 


13 CO 


2-1 


220.399 


12 




3.3 


e 


55 


0.2 


9.1 


4.3 x 10" 6 




5-4 


550.926 


14.7 




3.2 


c 


51 


0.5 


14.7 


1.1 x 10" 4 


cs 


2- 1 


97.981 


8.6 




3.0 


f 
J 


47 


0.2 


9.1 


3.5 x 10~ 7 




3-2 


146.969 


10.2 




3.0 


f 


20 


1.0 


10.0 


1.3 x lO" 6 




5-4 


244.936 


5.4 




3.0 


f 


13 


1.3 


5.5 


3.3 x 10- 6 




6-5 


293.912 


3.3 




2.2 


f 


12 


0.7 


3.1 


3.3 x 10~ 6 


[C I] 


1-0 


492.161 


3.7 




3.5 


A' 


54 


0.1 


3.7 


2.1 x 10~ 5 


[cm 


3/2-1/2 


1900.537 




4.4 x 10- 3 




h 


55 


1.0 


13.4 


4.4 x 10~ 3 


[O I] 


1-2 


4744.778 




4.1 x 10~ 3 




h 


46 


0.9 


0.88 


4.5 x 10" 3 



" All results for 7k = 55 K, n(H2) = 4x 10 s cm 3 . T\,= Rayleigh- Jeans brightness temperature, 4 nl = flux, Av = line wid th, T ex = excitation tem - 
perature, t = line-ce nt er peak opacity. b F rom H 2 0. c This work. d lUngerechts et all j!986l) . c iMinchin et"ai] d!993h . f ISnell et alj d!984al) . 
* iPlume et al.Ul994h . * lEmerv et alJ fl9961) . 



LTE and adding corrections for opacity and beam-filling, the to- 
tal source averaged column density can be calculated as 

C T Znkv-. 1 Q(T) f 
M.TE = — -ty 1 -r- exp (EJkT ex ) T b dv, (9) 
rjhi hc i A u i g u J 

where C t = t/(1 - exp(-r)) is the opacity correction factor. 
Using C T = 1, a source size of 85" to calculate tj\,{ and an excita- 
tion temperature of 10 K, the ammonia outflow column density is 
estimated to be ~1.2xl0 13 cirT 2 . This column density is thus not 
opacity-corrected and may therefore be underestimated. The am- 
monia outflow abundance is calculated as Z(NH3) = Nnh } /Nh 2 
and is found to be ~7xl(T 10 . 

4. Model results: RADEX - a homogeneous model 

In this section we use an enhanced version of the RADE30 
( Ivan der Tak et all [2007) code to compute the intensities of the 
narrow components of the transitions observed by Odin, and 
the intensities of several atomic and molecular transitions that 
previously have been observed in S140. The RADEX codex ap- 
plies a very simple method of mean escape probabilities to the 
radiative transfer and yields results that are similar to those in 
the large-velocity-gradient approximation, assuming an isother- 
mal and homogeneous medium, but without an implicit assump- 
tion of a gradient in velocity. Although much more sophisticated 
methods of radiative transfer have been used to construct models 
of the atomic and molecular line emission from this region, it is 
useful to examine a simple, internally consistent model of a ho- 
mogeneous cloud that reproduces the principal observed facts. 
Such a model can be used in particular to assess the role of 
radiative coupling of molecular excitation to the intense con- 
tinuous radiation within the PDR. In Sect. [5] these results are 

2 The published version will soon include these enhancements, see 
http : //www . sron . rug . nl/ ~ vdtak/r adex/index . shtml 



compared to the results from a 3D inhomogeneous PDR model 
that calculates the temperature and abundance distributions, and 
a multi-zone escape probability method that calculates the level 
populations and intensity distributions of water and ammonia 
JPoelman & Spaansll2005ll2006l) . 

The molecular data files for H2O and NH3 have been en- 
larged, compared with those previously available through the 
Leiden Atomic and Molecular Database (LAMDAfl Data for 
the ortho and para forms of these molecules have been com- 
bined into a single file and in frared transi t ions h ave been added. 
Collisi on rates published by iFaure et all d2007l) for H 2 + H 2 
and bv lFaure et al.l (2004) for e~ + H2O are now used. Two dif- 
ferent data files have been used to analyse the excitation of NH3. 
The smaller includes only the lowest 112 levels of the vibra- 
tional ground state and 140 radiative transitions between them. 
A larger file, derived mainly from the HITRA>Q database, has 
also been tested: it contains 2 392 levels and 15 067 radiative 
transitions, including lines of several v 2 and V4 vibrational bands 
in the mid-infrared. 

We include a simple dust model of t h e broa d -band spectrum 
dThronson et al.L 119831; IMinchin et all 119951; iNev & Merrill 
1980) at submm and far-infrared wavelengths in order to charac- 
terize the internal radiation sensed by the molecules (a detailed 
description is found in Appendix|A]of the on-line material). 

The principal parameters needed to specify a model 
are the kinetic temperature and the average number 
density of molecular hydrogen n(H2). Our adopted dust 
model implies /V(H2) =4.7xl0 22 crrT 2 and an average density 
n(H2) — 3xl0 4 cm" 3 in a homogeneous, spherical cloud with a 
source size L = 0.5pc (derivation in Appendix lAl on-line ma- 
terial). Unfortunately, this average density is inconsistent with 
our observations. Although a uniform RADEX model can be con- 
structed based upon this density, this low-density model cannot 

3 http : //www. strw. leidenuniv.nl/~moldata 

4 http : // cf a-www . harvard . edu/hi tr an 
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Table 3. Results of the homogeneous RADEX model for the nar- 
row PDR component: column densities and abundances." 



Species Column Density 


Abundance 


x N(x) 


X(x) 


[cm- 2 ] 





H 2 


2.0 x 10 14 


9.8 x 10~ 9 


H' 8 


< 4.9 x 10 12 


< 2.5 x 10" 10 


NH 3 


1.5 x 10 14 


7.5 x 10" 9 


I2 CO 


2.0 x 10 18 


1.0 x 10~ 4 


13 CO 


2.9 x 10 16 


1.4 x 10~ 6 


CS 


1.5 x 10 14 


7.5 x 10 -9 


C 


2.0 x 10 17 


1.0 x 10~ 5 


c + 


9.0 x 10 17 


4.5xl0~ 5 


o 


6.3 x 10 17 


3.2 x 10- 5 



All results for T K = 55K, «(H 2 ) = 4x 10 5 cm 4 , n(eT) = 7.7 cm" 3 , 

MH 2 ): 



sion transitions near 23 GHz. Indeed, the calculations indicate 
that the (J,K) = (3,3) inversion transition is a weak maser (cf. 
Walmsley & Unger echtslll983b . In addition, NH3 has strong vi- 
brational transitions in the mid- and near-infrared. When these 
are included, the intensity of the 572 GHz transition is somewhat 
suppressed while the lowest inversion transitions are somewhat 
enhanced. 



Table 4. Clumpy PDR model parameters. 



explain the observed intensities of the pure rotational lines of 
H2O and NH3 unless the line-center optical depths are of the 
order of 200 and 100, respectively. Such large opacities would 
imply significant line broadening through saturation of the emis- 
sion, which conflicts with the observed narrow profiles. 

High densiti es of 7xl0 5 cm" 3 a re found by multitransition 
CS observations dSnell et al.Lll984al). and 5xl0 5 cm 3 from NH 3 
inversion transitions (Ungerechts et all 1 1986b . We include the 
CS transitions in our RADEX model in addition to our own ob- 
servations, and our best fit suggests a mean molecular density 
of 4xl0 5 ctrr 3 and a temperature of 55 K. This suggests that the 
strongly emitting molecules occupy only a 7% fraction of the 
volume. The H2 column density, which is needed to obtain the 
abundances, X(x) = N(x)/N(U 2 ), is calculated from the jV( 13 CO) 
obtained with RADEX. This column density is also in agreement 
with AT ROT ( 13 CO) found in Sect. ED Using [ 12 C/ 13 C] =70 and 
[CO/H 2 ] = 10~ 4 we obtain N(U2) = 2.0xl0 22 ctrT 2 for the nar- 
row PDR component. 

Table|2]collects the observed intensities and the correspond- 
ing values calculated from the model. A source size of 120", 
derived in Sect. [2] is used to correct the antenna temperatures 
T* K to Rayleigh-Jeans brightness temperatures T\, according to 
the scaling in Eq. [5] The estimated abundances and column 
densities are summarized in Table [3] Note that the best-fitting 
model of the [C II] /I157yum line emission implies a density 
of carbon ions, n(C + ) = 7.7ctrT 3 . This is taken to be equal to 
the electron density, which implies an average electron fraction 
of X(e~) = 1.9xl0~ 5 . This electron fraction is large enough that 
electron-impact on polar molecules like H2O must be competi- 
tive with H2 collisions in exciting rotational states. However, the 
results presented for NH3 so far neglect electron collisions. 

The model reproduces nicely all main beam temperatures ex- 
cept the ammonia inversion transitions which are predicted to be 
lower than the observed intensities. The 13 CO(2 - 1) transition is 
also predicted to be slightly lower than what is observed. These 
discrepancies may be caused by a co-existing, extended lower 
excitation gas component, which is not contributing much to the 
emission of the 13 CO(5-4) or NH3(1q-Oo) lines. Even so, the 
colum n de nsities of C, 1 3 CO an d CS agree well with Plu me et alj 
d 19941) and lSnell et all dl984al). and the C abundance is also in 
agreement with Franklin et al J (20081). The w ater abundance was 
pr eviously estimated bv ISnell et alj d2b00a) to be 9xl0" 9 , and 
bv lAshbv et al] J2000a) to be 2xl0 -8 , both very similar to the 
abundance calculated by RADEX. 

The RADEX model is also tested for higher transitions of H2O 
and shows that some lines are sensitive to various effects. For 
example, the intensity ratio of the 1661 and 1670 GHz lines ap- 
pears to be sensitive to the total column density. When the inter- 
nal radiation field is switched off in the models, the intensities 
of the transitions involving the lowest rotational levels increase 
slightly while the intensities of lines involving more highly ex- 
cited states are diminished: in particular, the 1661GHz line is 
enhanced by radiative excitation. 

A series of excitation calculations with RADEX at small op- 
tical depths indicate that the intensity of the 22,1-2^2 line at 
1661GHz is very sensitive to the total column density: evi- 
dently photon-trapping is ineffective at column densities around 

10 14 ctrT 2 or less, so that this line would appear in absorption if ity dispersion of 1 .2 km s _1 dZhou et 
the opacity in the 557 GHz transition is small. 

The excitation of NH3 also shows some interesting effects 
in these models. Due to the rather high density of C + ions, the 
reaction of C + with NH3 will destroy ammonia at a rate ex- 
ceeding 10~ 8 s~', which is rapid enough to affect the excita- 
tion of the metastable levels involved in several of the inver- 



size 


"c(H) 


meCH) 


F" 




T 

1 gas 


^dust 


[pc] 


[cm- 3 ] 


[cm" 3 ] 


[%] 


[pc] 


[K] 


[K] 


0.5 


2xl0 6 


lxlO 4 


8 


0.03 


9-40 c 


18- 32 c 












35-195 rf 





"Volume filling fraction of clumps. 
''Clump size. 
'Temperatures in clumps. 
rf Temperatures in interclump medium. 



5. Model results: clumpy PDR model 

A self-consistent physi-chemical 3D inhomogeneous model is 
used to compute the temperature of the gas and dust, and the 
abundance distributions inside the clumpy cloud (detailed de- 
scriptions are found in ISpaansl [1996: Spa ans & van Dishoeckl 
1 19971: iPoelman & Spaansl 120051) . This clumpy PDR model has 
three free parameters; the volume filling factor F (fraction of 
the cloud that is occupied by clumps), a clump size l c , and the 
clump-interclump ratio n c (H)/«i c (H). The total hydrogen den- 
sity «h = «(H) + 2n(H2). In the dense clumps, hydrogen is en- 
tirely molecular, thus n(H) = 2n(H2). The PDR model parame- 
ters used in this paper are given in Table [4] W e use a veloc- 

al.L 119941) equivalent to a 
FWHM line width of 2.0 km s" 1 for an optically thin line. With 
an opacity of about 5-10 the line is broadened by a factor of 
~ 1 .5 — 2 to 3-4kms _1 , consistent with our observations. The 
inc ident radiation field is taken to be /uv = 140 with respect to 
the iDrainel (1 1 9781) field. The clumps are randomly distributed in 
the interclump medium and the volume filling fraction is taken 
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to be 8% dSpaans & van Dishoeckl [l997t lAshbv eTafl l2000d: 
lPoelman"& Spaans, 120051) . The grid size is 81x81x81, corre- 
sponding to a resolution of 0.006 pc and a total size of 0.5 pc 
(about 1 10" at a distance of 910pc). 

In the PDR computation, each clump is joined smoothly to 
the inter-clump medium through a power-law ocl/R 2 density dis- 
tribution. This envelope contains very little molecular material 
and the bulk of the impinging radiation field is always absorbed 
by the clumps and not the inter-clump medium. The definition 
of the mean H2 density and the mean molecular abundances X 
of species x is thus made through < n >- Fn c + (\ - F)n lc 
and X mean (x) = F X c (x) + ( 1 - F) X ic (x) to allow an unambiguous 
comparison with molecular line observations . 

A high molecular hydrogen density of ^lxlO 6 cm 4 in the 
clumps is required to match the intensities observed by Odin 
and SWAS. Note th at this clump densi t y is consider ably higher 
than that used by iPoelman & Spaansl d2005l 120061) . The high 
density and clumpy medium suggested by our PDR model is, 
however, also supported by previous density estimates. For ex- 
ample, density variations between 10 4 - 10 6 ctrT 3 are necessary 
to explain the apparent coincidence of the CI and CO ridges 
as well as the distance betw een the CI ridge and the bright 
rim (Hayashi&Murata, 1992). Emission from numerous, small, 
dense (10 5 -2xl0 6 cm 4 ) and optically thick clumps of gas, 
where the number density of clumps decreases with distance 
from the CS cloud centre close to IRS, is the best explanation of 
the CS column density variations, as well as the observed opaci- 
ties an d temperatures obta i ned from multitra nsition observations 
of CS (ISnell et all Il984ab IZhou et all 11994 . As molecular hy- 
drogen density in the interclump medium we use 5xl0 3 cm 4 . 



log(intensity 0- H 2 O) 



erg 



s cm' sr 



120 



100 



80 



60 



40 



20 




I 

20 



) 60 
offset [arcsec] 



~~r~ 

80 



-10 



-12 



100 



Fig. 6. Result of the PDR high-density model. The water inten- 
sity is shown in a logarithmic scale. The ionizing star is located 
on the left side of the displayed region. 

Table 6. Observed and predicted T* A in the PDR for 
o-H20(li o- 1q,i) with SWAS and Odin at the central position. 



Table 5. Resulting abundances with respect to «h 2 for the narrow 
PDR components of water and ammonia from the clumpy PDR 
model. 



Species 


Comp. 


m 2 


Abundance 






[cm" 3 ] 




H 2 


Interclump 


5.0xl0 3 


8.9x10-'° 




Clump 


l.OxlO 6 


4.8xl(T 8 




Mean 


8.5xl0 4 


4.6x10-' 


NH 3 


Interclump 


5.0xl0 3 


1.1x10-'° 




Clump 


l.OxlO 6 


4.5xl0- 8 




Mean 


8.5xl0 4 


3.7xl0- 9 



The resulting abundances, with respect to H2, in the clumps 
and the inter-clump medium, as well as a mean value over the 
whole region, are found in Table [5] The mean value is about 
half than found by RADEX. In the clumps, the water abundances 
vary between lO -7 - lO -9 with an average of 4.8xl0- g . The low 
density in the interclump medium results in low water abun- 
dances between 1x10-" -5xl0- 9 with an average of 8.9xl0 10 . 
Ammonia is found to have about the same average abundance as 
water, 4.3 xlO 8 , in the clumps, while lower by almost an order 
of magnitude in the inter-clump medium, 1.1x10 . 

The result from our PDR model is used as input to a three- 
dimensional multi-zone escape probability method which com- 
putes the level populations of ortho- and para-fyO (up to 
~35 0K), CO and NH3, as well as line intensities and opaci- 
ties (IPoelman & Spaansl 2005, 2006). The excitation due to dust 
emission is fully included. The dust temperature in the model 
varies between 18 and 32 K. However, due to geometry the high- 
est Tdust is only achieved at the edge of the PDR and the dust 



Satellite 


obs. r; 


Pred. T* 




[mK] 


[mK] 


SWAS 


186 


120 


Odin 


416 


345 



component fills only a small part of the region. The bulk of the 
dust has a temperature of about 25 K and the dust continuum 
opacity /UOOyum is 0.02 or smaller. Results from this code dif- 
fer at most by 10% compared to other Monte Carlo/Accelerated 
Lambda Iteration codes. Fig. [6] shows the predicted water inten- 
sities on a logarithmic scale. Note that the ionizing star is located 
on the left side of the displayed region. The peak emission is 
found about 70-80" from the bright rim, in agreement with our 
observations of all species. 

The resulting model intensities are then used to calculate the 
expected antenna temperature via 



1 SAeff 

2 kAv 
1 bA 2 



1 br] b{ Q mb A 2 ^ 

2 kAv D. A 



2 kAv 



?7bf T]mb 



kAv 

bA 3 

2*A^ mb%f [K] ' 



(10) 



where S is the integrated flux density in [erg s cm ], the 
intensity b is in [ergs -1 crrr 2 sr -1 ], A e ff is the effective antenna 
area 77A A geo = A 2 /Q.a [cm 2 ], Av is the observed line width in Hz, 
fi m b is the main beam solid angle in [sr], Av is the observed 
line width in cms 1 , and the wavelength A is in cm. We con- 
volve the intensities to the Odin beam, and in addition predict 
the intensities of the forthcoming observations of the Herschel 
Space Observatory using a telescope diameter of 350 cm and a 
beam-efficiency of 0.7. The observed and modelled SWAS and 
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Fig. 7. Predicted antenna temperatures of the ortho-WiOilx $ — 
lrj,i) and para-HzOil ij -Oo,o) transitions for the forthcoming 
Herschel-HIFI observations along a cut 20-100" from the 
bright rim. 



Odin antenna temperatures for the l^rj- lo,i transition are pro- 
vided in Table [6] for the central position. In addition, the pre- 
dicted antenna temperatures for Herschel-HIFI of both the or- 
tho- and /5ara-H20 li,rj — lo,t an d li.i — Oo,o transitions along a 
cut from the bright rim are shown in Fig. [7] 

The model ortho-water opacities are relatively low with a 
unweighted mean opacity of 17 over the entire source, with a 
range between 10~ 5 - 800. 

The antenna temperatures from the PDR model agree well 
with our observations, especially considering uncertainties in 
source and clump size, as well as the volume filling fraction, 
which is the main uncertainty and can be varied to obtain slightly 
higher or lower T* K . The dust background continuum is also 
lower than in the RADEX model, which diminishes the effect 
of radiative excitation of water. A clumpy medium with high 
clump densities above 10 6 ctrT 3 is, however, necessary to pro- 
duce model antenna temperatures close to the observed val- 
ues. This is also supported by the resulting low density of 
3xl0 4 cm~ 3 obtained by the dust model in Sect. [4] assuming a 
homogeneous, spherical cloud, which is about 7% of the density 
in our best fitting model. 



6. Summary and conclusions 

We have used the Odin satellite to observe water, ammonia and 
carbon-monoxide in the well-known molecular cloud S140. We 
have simultaneously observed five-point strip maps across the 
bright rim in S140 of the ortho-^O li,o-lo,i and the ortho- 
NH 3 ln-Oo transitions, as well as CO(5-4) and 13 CO(5-4). 
The NH3 transition has never previously been observed in S140. 
Observations of Hi O in the central position resulted in a non- 
detection at a rms level of 8 mK. As support observation we also 
mapped 13 CO(l -0) with the Onsala 20m telescope. 

Like CO, the water line-profile is dominated by emission 
from a NW-SE outflow, however, mainly in the red wing. 
Strong self-absorption is seen in the optically thick CO emission, 
while no obvious signs are seen in the ortho-fyO, ortho-Nil?, or 
the almost optically thin 13 CO line profiles. In addition to the 
outflow, our water line shows emission from a more extended 
NE - S W elongated PDR. Both these components originate ap- 
proximately around our central position. No additional emission 
closer to the bright rim or further into the molecular cloud is 



detected. Close to the bright rim the temperature is most likely 
too high for a detection of our transition with an upper state en- 
ergy of 61 K. Instead, higher-lying transitions, observable with 
the Herschel Space Observato ry, will have their peak intensity 
shifted towards the bright rim (Poelman & Spaani|2005j,|2()06). 
Even closer to the bright rim at a few magnitudes of Ay, water 
is, however, photo-dissociated by the UV field. 

The ortho-NHy, emission seems to emanate from the same 
high density clumps in the PDR and the outflow as water, but 
also shows additional emission further into the cloud where the 
ambient gas temperature drops to about 30 K. Compared to water 
in the central position, ammonia has a weaker outflow emission 
in the red wing although similar in the blue outflow. Close to the 
bright rim, where the outflow contribution to the emission is very 
low, the water and ammonia line profiles are almost identical, 
suggesting an origin in the same gas and velocity fields of the 
PDR. The 13 CO line also shows a very similar line profile as H2O 
and NH3 in this position with a narrow line width of ~3 km s _1 . 

Abundances, with respect to H2, in the PDR are estimated 
both with an enhanced version of the homogeneous RADEX code 
and with a clumpy PDR model. This model points to low mean 
water and ammonia abundances in the PDR of 5xl0~ 9 and 
4x1 0~ 9 , respectively. In the high-density clumps both the av- 
erage water and ammonia abundances increase to 5xl0^ 8 . To 
match the observed PDR antenna temperatures with Odin and 
SWAS, a clumpy medium is required by the model with a high 
molecular hydrogen density in the clumps of ^ 1 x 1 6 cm -3 . The 
resulting RADEX mean abundances are twice as high: 1.0xl0~ 8 
and 8xl0~ 8 for water and ammonia, respectively, using a molec- 
ular hydrogen density of 4xl0 5 cirT 3 and a kinetic temperature 
of 55 K. The differences most likely arise from the uncertainty 
in density, beam-filling, and volume filling of the clumps. The 
opacity of the narrow PDR component of the 1 10 — lo,i transi- 
tion is constrained by the narrow line width, and is estimated by 
RADEX to be ~7. The PDR model also confirm a low water opac- 
ity with an unweighted mean opacity of 17 and a model range of 
~10" 5 -800. The mean outflow water abundance, derived from 
a simple two-level approximation, is higher than in the PDR by 
at least one order of mag nitude, ~2xl0~ 8 -2xl0 -7 . 

Predictions of antenna temperatures for observations with 
HIFI are given by our PDR model of the ortho- and para-HiO 
1 1,0 - 1q,i and li,i — Oo ; o transitions for nine positions across the 
bright rim, and are found to peak around 70 - 80" from the dis- 
sociation front in agreement with our observations. 
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Appendix A: RADEX - construction of a dust model Appendix B: Figures and Tables 



In order to relate the molecular column densities, N(x) of species 
x, to fractional abundances, X(x) — N(x)/N(H2), a uniform, ho- 
mogeneous sphere of diameter L = N(H2)/n(H2) is assumed 
here. The adopted physical diameter of the PDR corresponding 
to an angular diameter of 120" (Sect. [3]) at a distance of 910pc 
is 0.53 pc. This is assumed to be equal to the line-of-sight depth. 

The observed intensity of the continuum is used to estimate 
the internal radiation field sensed by the molecules. We construct 
a simple model of the broad-band spectrum at submm and far- 
infrared wavelengths in order both to characterize the internal ra- 
diation and to estimate the tota l column densities of dust and hy- 
drogen. iThronson et a il (119831) measured the far-infrared emis- 
sion of S 140 and found a peak flux density of the order of 10 4 Jy 
slightly shortward of /H00/./m in a 49" beam. Minchin et al. 
( 1 1995b presented total broad-band fluxes in a 1.'5 x 1.'5 box. We 
represent the latter results with a two-component model of ther- 
mal emission by dust over a solid angle of Q. = 1.9 x 10~ 7 sr. 
The main component has a dust temperature 7d ust = 40 K and a 
long-wavelength (A > 40 /mi) form of the opacity law 

Tdust = 0.0679(100M) L2 

where A is the wavelength in /urn. The opacity law is smoothly 
matched to a standard interstellar extinction law at shorter wave- 
lengths, which is also used to describe the second component 
at Tdust = 140 K. The opacity of the first dust component corre- 
sponds to a visual extinction A y = 58.8 mag. The second compo- 
nent has a smaller optical depth Ay = 0.023 mag, but is assumed 
to cover the s ame solid angle . In ad dition, the mid-infrared mea- 
surements of Ney & Merrilll Jl980h have been adapted in order 
to specify the radiation field at even shorter wavelengths. In the 
calculations, the molecules are assumed to be exposed to an av- 
erage intensity of continuous radiation 



b v (t C mb) + n 



J V 



n 



where B v is the Planck function, Tcmb = 2.73 K is the tempera- 
ture of the cosmic background radiation, / v dust is the flux density 
of the 2-component dust model, Q. = 1.9xl0~ 7 sr, and r\ = 0.72 is 
a dilution factor to scale the brightness of the dust source to the 
larger beam area of the Odin measurements. It is important to 
keep in mind that we observe this strong far-infrared radiation; 
therefore, the co-extensive molecules must sense it also. 

For the adopted interstellar extinction law and a standard 
gas/extinction ratio, 

2/V(H 2 ) = 1-6 x 10 21 A V cm" 2 , 

the adopted dust model implies NQI2) =4.7xl0 22 ctrT 2 and 
an average density n(H2) = 2.9xl0 4 cirT 3 over the source size 
L = 0.53pc. This average density is, however, inconsistent with 
the observed molecular line emission in large beams (9 > V). 
Although a uniform RADEX model can be constructed based 
upon this density, the line-center optical depths of the pure 
rotational lines of H2O and NH3 would be of the order of 
200 and 100, respectively. Such large opacities would imply 
significant line broadening through saturation of the emission, 
which conflicts with the observed narrow profiles of ~3 km s~' . 



547.676 GHz 




-30 -20 -10 

Velocity [km/s] 



Fig. B.l. Odin observations of Hl s O in the central position. 



6 Broad component 




Upper state energy [K] 



Fig. B.2. Rotation diagram of the broad component of 13 CO(l - 
0) w ith the Onsal a 20-m telescope, J = 2-1 and 7 = 3-2 
from Min chin et al.l d 1993b . and 7 = 5-4 with Odin, producing 
r R0T = 24+2 K and /Vrot = (2.5±0.4)x 10 16 cm" 2 . 
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Table B.l. Observed transitions and their parameters" in S140 with the Odin satellite in a five point NE-SW strip. 



Species 


Transition 


Freq 
[GHz] 


E u 
[K] 


A-coeff 
[s- 1 ] 


Av 

[kms- 1 ] 


Pos. 


U LSR 

[kms-'] 


[mK] 


T : 

J n& 

[Kkms- 1 ] 


rms 
[mK] 


note 


H 2 


1 1,0 - 1(U 


556.936 


61.0 


3.46e-3 


0.27 


1 


-8.0 


207 


0.64 


39 
















2 


-6.9 


517 


2.74 


68 
















3 


-6.6 


632 


3.34 


19 
















4 


-7.0 


379 


1.84 


57 
















5 


- 


- 


- 


44 


No detection. 


H' 8 


1 i,o — lo.i 


547.676 


60.5 


3.29e-3 


0.34 


3 


- 


- 


- 


~8 


No detection. 


CO 


5-4 


576.258 


83.0 


1.22e-5 


0.32 


1 


-5.4 


6 680 


47.1 


67 
















2 


-5.0 


11760 


94.2 


47 
















3 


-5.5 


13 240 


108.8 


43 
















4 


-6.6 


8 690 


60.1 


50 
















5 


-6.9 


2 220 


10.1 


160 




13 CO 


5-4 


550.926 


79.3 


1.10e-5 


0.27 


1 


-7.8 


2 740 


8.29 


115 
















2 


-7.4 


5 930 


21.4 


81 
















3 


-7.4 


7 110 


27.7 


61 
















4 


-7.6 


3 640 


13.2 


119 
















5 


-8.5 


1 190 


2.5 


291 




NH 3 


lo-Oo 


572.498 


27.5 


1.61e-3 


0.32 


1 


-8.0 


444 


1.18 


99 
















2 


-7.4 


467 


1.83 


33 
















3 


-7.4 


588 


2.66 


29 
















4 


-7.3 


381 


1.86 


36 
















5 








91 


No detection 



" Transition = the quantum numbers for the transition; Freq = rest frequency of the transition; E a = the upper state energy; Acoeff = the Einstein 
A-coefficient; Av = the velocity resolution; Pos = the strip position from NE to SW. v LSR = the peak LSR velocity; r pcak = the observed peak 
temperature of the transition; j T^dv = the integrated intensity from the observed spectra not corrected for beam-efficiency or beam-filling; 
rms = noise. 



33 3 




Upper state energy [K] 

Fig. B.3. Rotation diagram of the narrow comp onent of 
13 CO(2-l) and 7 = 3-2 from iMinchin et al.l d 19931) . 7 = 5- 
4 with Odin, and 7 = 6-5 from iGraf et all d 19931) . producing 
r ROT = 69+27K and N RO j = (3.2+1. 8)xl0 16 cm" 2 . 13 CO(1-0) 
is not included in the fit. 




-30 -30 -10 10 20 

Velocity [km/s] 

Fig.B.4. Gaussian fits to 13 CO(5-4) at the central position. 
The widths, amplitudes and centre velocities are 3.2kms _1 and 
8.2kms -1 ; 6.610K and 0.612K; -7.3 kms 1 and -6.8kms ', 
respectively. 
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20 



Velocity [km/s] 

Fig. B.5. Gaussian fits to the convolved (to the Odin 126" beam) 
spectra of 13 CO(1-0) at position 1. The widths, amplitudes 
and centre velocities are 2.7 km s -1 and 8.6kms~'; 6.982 K and 
0.395 K; -7.6kms~' and -8.0kms~', respectively. 




Velocity [km/s] 

Fig. B.6. Gaussian fits to H2O at the central position. The widths, 
amplitudes and centre velocities are 3.1kms~' and 8.8kms _I ; 
416 mK and 213 mK; -7.1 kms -1 and -6.1 kms -1 , respectively. 




-30 -20 -10 10 



Velocity [km/s] 

Fig. B.7. Gaussian fits to NH3 at the central position. The widths, 
amplitudes and centre velocities are 3.3 kms -1 and 8.5 kms -1 ; 
487 mK and 100 mK; -7.6 kms -1 and -6.4 kms -1 , respectively. 
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Table B.2. 13 CO Gaussian fits". T h uses a source size for the PDR (narrow component) of 120"— * rj^f = 2, and a source size for the 
broad outflow component of 85"— > rj^f = 3. 



Pos. 


fLSR 


err 


T A , 


T b 


err 


Av 


err 




[kms- 1 ] 


[kms" 1 ] 


[mK] 


[mK] 


[mK] 


[kms" 1 ] 


[km s" 1 ] 


1 


-7.7 


0.04 


2 665 


5 920 


68 


2.8 


0.08 


2 


-7.3 


0.02 


5 877 


13 060 


56 


3.3 


0.04 


3 


-7.3 


0.01 


6610 


14 690 


199 


3.2 


0.07 


3 


-6.8 


0.33 


612 


2040 


207 


8.2 


1.54 


4 


-7.3 


0.03 


3 620 


8 040 


66 


3.5 


0.07 


5 


-8.2 


0.25 


1 124 


2 500 


282 


2.1 


0.62 



" Gaussian fits (including errors) to the spectra in 5 positions. For most positions two Gaussians are needed to fit the spectra. Parameters: 
Vlsr = the LSR velocity at the peak temperature, amp = the amplitude of the Gaussian fit, Av = FWHM line width. 



Table B.3. H2O Gaussian fits . 7b(PDR) uses a source size of 120"— > tjm = 2, while 7b(outfiow) uses a source size of 85"— » t/m = 3. 



Pos. 


VLSR 


err 


T A , 


T b 


err 


Av 


err 




[kms" 1 ] 


[kms" 1 ] 


[mK] 


[mK] 


[mK] 


[kms" 1 ] 


[kms" 1 ] 


1 


-7.5 


0.10 


180 


400 


12 


3.1 


0.2 


2 


-6.9 


0.08 


340 


756 


28 


3.5 


0.3 




-5.4 


0.39 


171 


570 


27 


8.8 


0.7 


3 


-7.1 


0.06 


416 


924 


30 


3.1 


0.20 




-6.1 


0.24 


213 


710 


31 


8.8 


0.60 


4 


-7.1 


0.08 


316 


702 


44 


3.5 


0.4 




-6.6 


0.8 


78 


260 


45 


10.0 


4.1 



" Notation as in Table |R2l. 

Table B.4. NH3 Gaussian fits . 7b(PDR) uses a source size of 120"— > rj\,f = 2, while 7b(outfiow) uses a source size of 85"— > T]\,f = 3. 



Position 


^LSR 


err 


T A - 


T b 


err 


Av 


err 




[kms- 1 ] 


[kms" 1 ] 


[mK] 


[mK] 


[mK] 


[kms" 1 ] 


[kms- 1 ] 


1 


-7.5 


0.14 


438 


973 


46 


2.6 


0.3 


2 


-7.9 


0.25 


A\2 


916 


200 


3.3 


0.9 




-6.1 


3.4 


111 


370 


160 


5.8 


2.6 


3 


-7.6 


0.12 


487 


1080 


90 


3.3 


0.50 




-6.4 


1.4 


100 


333 


92 


8.5 


4.3 


4 


-7.4 


0.23 


321 


713 


185 


3.9 


1.2 




-6.8 


2.6 


61 


203 


192 


8.7 


12.1 



" Notation as in Table HOI. 



